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Abstract 

Recent experiments revealed that the dielectric dispersion spectrum of fission 
yeast cells in a suspension was mainly composed of two sub-dispersions. The 
low-frequency sub-dispersion depended on the cell length, whereas the high- 
frequency one was independent of it. The cell shape effect was qualitatively 
simulated by an ellipsoidal cell model. However, the comparison between 
theory and experiment was far from being satisfactory. In an attempt to close 
up the gap between theory and experiment, we considered the more realistic 
cells of spherocylinders, i.e., circular cylinders with two hemispherical caps at 
both ends. We have formulated a Green's function formalism for calculating 
the spectral representation of cells of finite length. The Green's function can 
be reduced because of the azimuthal symmetry of the cell. This simplification 
enables us to calculate the dispersion spectrum and hence access the effect of 
cell structure on the dielectric behavior of cell suspensions. 

PACS Numbers: 82.70.-y, 87.22.Bt, 77.22. Gm, 77.84.Nh 
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I. INTRODUCTION 



Dielectric spectroscopy has become a quantitative method of real-time monitoring of cell 
growth in suspensions [1-3]. The real-time monitoring has advantages over conventional 
techniques and would be important to investigate the dynamic behavior of cell growth. 
There are many factors that may influence the dielectric behavior of biological materials: 
structure, orientation of dipoles, surface conductances, membrane transport processes, etc. 
All these factors influence one another and it is difficult to separate out the effect of a single 
one. However, some effects can be dominant at certain ranges of frequencies. For instance, 
the dielectric dispersion spectrum of living cell suspensions in radiofrequencies is mainly 
determined by the cell shape. The objective of this work is to develop a theory for such 
correlation, on which new applications in biotechnology rely. 

More recently, Asami reported that the dielectric dispersion spectrum of fission yeast 
cells in a suspension was mainly composed of two sub-dispersions. The experimental data 
revealed that the low-frequency sub-dispersion depended on the cell length, while the high- 
frequency one was independent of it. Asami adopted a shell-ellipsoid model [3], in which an 
ellipsoid is covered with an insulating shell as the electrical model of nonspherical biological 
cells. The comparison between model calculation [3] and experimental data |^ was far from 
being satisfactory. Asami suggested that the discrepancy is attributed to the rod-like cell 
shape. However, the depolarization factor needed in his theory is difficult to calculate for 
cells of rod-like shape because of lack of available theories. 

In this work, we propose the use of the spectral representation P| for analyzing the cell 
models. The spectral representation is a rigorous mathematical formalism of the effective 
dielectric constant of a two-phase composite material 0. It offers the advantage of the 
separation of material parameters (namely the dielectric constant and conductivity) from the 
cell structure information, thus simplifying the study. From the spectral representation, one 
can readily derive the dielectric dispersion spectrum, with the dispersion strength as well as 
the characteristic frequency being explicitly expressed in terms of the structure parameters 
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and the materials parameters of the cell suspension (see section III below). The actual 
shape of the real and imaginary parts of the permittivity over the relaxation region can be 
uniquely determined by the Debye relaxation spectrum, parametrized by the characteristic 
frequencies and the dispersion strengths. So, we can study the impact of these parameters 
on the dispersion spectrum directly. 



II. FORMALISM 

A. Spectral representation theory 

We consider a two-component composite dielectric with complex dielectric constant e(r) 
equal to e2 = £2 + o"2/j27r/ in the host medium and ti = £1 + (Ji/j2nf in the embedding 
medium. An interface S separates the two media. In a uniform applied electric field E = Eqz 
(for convenience, let Eq = —1), the electrostatic potential satisfies the Laplace's equation: 

V ■ [(1 - %(r))V$(r)] = 0, (1) 

where s = 62/(62 — ei) denotes the relevant material parameter and ri{r) is the characteristic 
function of the composite, having value 1 for r in the embedding medium and otherwise. 
The electric potential $(r) can be solved formally as: 

$(r) = <l>o(r) + dr'r]{r')VGo{r - r') ■ V'<l>(r') (2) 

where Go{r — r') = l/47r|r — r'| is the free space Green's function, and $o(r) = r • z = z 
is the potential of the unperturbed uniform field E. It is instructive to convert the volume 
integration into the surface integration via the Green's second identity and only deal with 
the potential on the interface S 0. We denote an integral-differential operator F: 

r<l>(r) = / dS' ■ V'Go(r - r')<I>(r') + -<^{r), r G S, (3) 
Js 2 

to avoid the singularity of G'o(r— r') when the integration variable r' approaches the point of r 
1^. The integration with a "prime" denotes the restricted integration which excludes r' = r. 
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Let \I/n(r) and s„ be the nth eigenfunction and eigenvalue of the F operator respectively. We 
can expand $o(r) and ^(r) on the interface S in a series expansion of eigenfunction \l/„(r): 



$o(r) = ^2;„^„(r), 

n 

$(r) = E^^n(r) 



(4) 
(5) 



where z„ are the expansion coefficients. Then we can write the effective dielectric constant 
e in the Bergman-Milton representation 

e = -l| dVe{v)E-, 

,(9$ 



I r 1 (9$ 

_/,V.,|l--,(r)]^ 



£2 



£2 



1--E 



dS ■ z$Jrl 



fn 



S - Sr 



(6) 
(7) 
(8) 

(9) 



where p is the volume fraction of the suspending cells. Note that Ez is a dimensionless 
electric field because Eq = —1. The eigenvalue s„ and the spectral function /„ can be 
proved to be real and satisfy simple properties that < s„ < 1, /„ > and J2 fn = 1- We 
shall show that both the spectral function and the eigenvalue s„ determine the dielectric 
behavior of cell suspensions. 



B. Cells with an axis of revolution 

Now the principal problem is to calculate the eigenvalue s„ and the spectral function /„. 
For many cells interacting with one another, it is a formidable task. However, in the limit 
of a dilute cell suspension and weak applied field, one can regard the cells in suspension 
as being noninteracting and randomly oriented and the problem is reduced to that of a 
single cell. We will consider cells with an axis of revolution, namely, the spheroidal and the 
spherocylinder cells [7] to mimic cells of rod-like shape. The prolate spheroid is generated 
by rotating an ellipse around its major axis, while the spherocylinder is obtained by fitting 
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two hemispherical caps at both ends of a circular cylinder. For a prolate spheroid, the 
eigenvalues and eigenfunctions can be calculated exactly. The only nonzero equals unity 
for E being along the major or minor axis of the prolate spheroid, and the corresponding 
eigenvalues are given by: 

1 +iii%±(4^, (10) 



' (g2_ 1)3/2 

= (1 - s,)l2. (11) 

where q is the length to diameter ratio; z and x refer to the direction along the major and 
minor axis respectively. For a spherocylinder, the consideration of the symmetry properties 
of the cell will help us choose the appropriate orthogonal basis for calculating the matrix 
elements of the T operator. Because of the rotation symmetry about the major axis of the 
spherocylinder, the eigenfunction is necessarily of the form (a„cosn^ + bnSmn9)f{x) with 
n being an integer. Due to the inversion symmetry of the cell, f{x) must be either odd or 
even functions. It is convenient to expand f{x) in a series of Legendre polynomials Pm{x/l), 
where 21 is the length of the cell. The applied uniform field E can always be resolved into two 
components along the major and minor axes of the spherocylinder, so we can calculate the 
Sn and /„ for E along the major and minor axes separately. By symmetry, in order to obtain 
a nonzero /„, the eigenfunction should be the form of ^ AjnP2m+i{x / 1) for E being along 
the major axis, while it reads cos6^BmP2m{xll) for E being along the minor axis, with 
m = 0, 1, 2, • • •. Using this orthogonal basis, we can calculate a truncated matrix according 
to the precision needed. We should remark that the matrix is generally nonsymmetric. 

III. DIELECTRIC DISPERSION SPECTRUM 

We show here that from the spectral representation, one can readily derive the dielectric 
dispersion spectrum. Substituting ei = £i + (Ti/j27r/ and 62 = £2 + C2/j27r/ {e and a being 
the real and imaginary parts of the complex dielectric constant) into Eq.(9), defining a new 
parameter t = (J2/{o'2 — o'i) and re-defining s = £2/(^^2 — ^1), we rewrite the effective dielectric 
constant e after simple manipulations: 



where en and are the high-frequency dielectric constant and the low-frequency con- 
ductivity respectively, while Ae„ are the dispersion magnitudes, are the characteristic 
frequencies of the nth sub-dispersion. 

We have already shown that there are only two poles in the spectral representation of 
the prolate spheroids. In what follows, we will show that there are two dominant poles in 
the spectral representation of the spherocylinder and hence there are two sub-dispersions in 
the dielectric dispersion spectrum. The dispersion magnitudes and dispersion frequencies 
are given by: 

3 s[s — si)[t — siY 
2 S2{t-sf 

- 2ne2tis - s^y ^^^^ 

Thus, we are able to obtain the dispersion strengths as well as the characteristic frequencies 
explicitly in terms of the structure parameters and the materials parameters of the cell 
suspension. 

The dielectric dispersion spectrum of a dilute suspension of prolate spheroids is mainly 
composed of two sub-dispersions, namely, Sz is responsible for the lower frequency one and 
Sx for the higher one. For a spherocylinder, we obtain a nonvanishing series of /„ and s„. 
Along the major axis, /i is dominant for all q and we can omit the smaller ones. This 
dominant /i is plotted in Fig.|I] against q, and the corresponding are plotted in Fig.|^, 
together with the exact result of a prolate spheroid. As is evident in Fig.[l|, we can see that 
the difference between the two models is indeed small. Along the minor axis, the solution 
becomes more complicated. The dominant /2 near q = 1 decreases quickly as q increases; 
another increases and takes over at large q. These two /„ are also plotted in Fig.|l| and 
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their corresponding eigenvalues are plotted in FigJ^. As shown in Fig.|^, the two eigenvalues 
tend to that of a prolate spheroid in the limit of both small and large q. 

Near q = 2, the two /„ have comparable values, resulting in two sub-dispersions at higher 
frequency. These sub-dispersions can interfere with each other, rendering it difficult to find 
the characteristic frequencies of the different sub- dispersions. Physically, the local field is 
the most nonuniform in this case. Nevertheless, we will consider cells of large length and 
omit this complication. 

With Eqs.(13)-(16), it is easy to calculate the effect of the rod-like cell structure on the 
dispersion spectrum and to compare with experiment data. We will show that the sphero- 
cylinder model does give some improvement towards the experimental result. However, the 
improvement is too small to close up the gap as Asami expected. In fact, we shall see that 
Asami omitted the material parameters which will play an important role in the experimen- 
tal condition. By introducing the conductivity contrast t = (T2/(o"2 — we found that a 
small negative t, i.e. o"i ^ 0-2, should be used to close up the discrepancy. 

We estimate t and s by fitting Eqs.(13)-(16) to the experimental ratio of Aei/Ae2 and 
fi/fu S^^ t — —0.0014 and s = 5.0. It means that cti ~ 700ct2 and ei ~ 0.80e2- 

The enhanced conductivity of cell cytoplasm is attributed to the membrane potential. The 
result is in contrast to the previous (unjustified) claim that ai ~ a"2. 

Table | lists the Aei/Ae2 ratio and /2//1 ratio for both experimental and theoretical 
results. Using the fitting material parameters, the improvement is obvious for both the 
prolate spheroid model and the spherocylinder model, while the difference between the two 
models is quite small. 

DISCUSSION AND CONCLUSION 

In this work, we have applied the spectral representation to the dielectric dispersion of 
suspensions of fission yeast cells. As mentioned by Asami the discrepancies between 
theory and experiment may be attributed to the rod-like cell shape. For cells of noncon- 
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ventional shape, however, there exists no available cell model in the literature and we must 
develop the spectral representation from first principles. 

More precisely, we have developed a Green's function formalism [6,8] for calculating 
the spectral representation of rods of finite length. We modelled the rod-like cells as the 
spherocylinders, i.e., circular cylinders with two hemispherical caps at both ends. We solved 
the spectral representation of the effective dielectric constant from first principles. Similar 
formalism was adopted for cell suspensions near their sub-division point [9-11]. 

Generally speaking, when the axial ratio q is larger than 4, the prolate spheroid model can 
be employed as a good approximation for rod-like cell structures. For g < 4, the dielectric 
behavior will become more sensitive to the cell structure of the suspending particles, and 
there are more sub-dispersions than that of the prolate spheroid suspensions. 

Our model does not include the rotational or vibrational effects, and our results are 
expected to be valid only for weak electric fields. Our model also ignores the multi-shell 
nature of the cells. Usually the multi-shell model is used to explain the high-frequency 
steps of spherical cell suspensions. Similar conclusions were found in one of our previous 
paper on multi-shell dielectric spheres in electrorheological (ER) fluids []T2[ to account for 



the effects of water coating on the ER effects. In Ref. WM, we also showed that the spectral 



representation can still be used for multi-shell model, albeit with a slight modification. 

In Asami's experiment, there exist three sub dispersions, the highest frequency step 
(above lOMHz) is due to the vacuole and cell wall as mentioned by Asami, while the two 
lower frequency step is evidently dependent on the cell shape. And the disperson magnitude 
of the highest frequency step is much smaller than that of the two lower frequency ones. 
So one expects that the multi-shell model has only a small effect on the lower frequency 
steps. In fact the multi-shell model was used in Asami's theory, but the discrepancy, as we 
mentioned in the text, is still large. 

The large cytoplasmic conductivity is a key result of our investigation. We believe that 
the large cytoplasmic conductivity is reasonable because the cells have to maintain a higher 
ion concentration in their cytoplasm to avoid the shrinkage of cells due to a loss of water 
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across the cell membrane. However, to our knowledge, there exists no direct experimental 
mesurement on the cytoplasmic conductivity. In our work, we propose a convenient and 
practical means of determining the cytoplasmic conductivity from the dielectric spectroscopy 
data, which analysis can be important for biotechnology. 
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TABLES 

TABLE I. The ratios of the characteristic frequencies Aei/Ae2 and the ratios of the dispersion 
strengths /2//1 hsted as a function of the length to diameter ratio q of the cells. The exper- 
imental results were extracted from Ref.[4] together with the theoretical predictions. Both the 
prolate spheroid model and the spherocylinder model adopt the same fitting material parameters 
determined from the experimental data. 

prolate spheroid spherocylinder 
experimental result Asami theory model model 

q Aei/Ae2 flj/fj Aei/Aea flj/fj Aei/Ae2 flj/fj Aei/Aea fj/R 
3.46 2.22 8.95 0.900 3.00 2.26 5.34 2.77 5.89 

7.17 8.65 27.4 2.07 7.73 6.10 15.3 6.67 16.4 

10.2 16.4 52.6 3.39 13.0 9.94 25.9 10.5 27.2 
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FIGURES 

FIG. 1. The three dominant /„ plotted against the axial ratio q: f\ along the major axis (solid 
line), /i along the minor axis (long dashed line) and /2 along the minor axis (short dashed line). 

FIG. 2. The eigenvalue s„ associated with along the major axis plotted against the axial 
ratio q: the spherocylinder cell (solid line with filled circles), and the exact result of the prolate 
spheroid (solid line). 

FIG. 3. The eigenvalue s„ associated with /„ along the minor axis plotted against the axial 
ratio q: si of the spherocylinder cell (solid line with filled squares), S2 of the spherocylinder cell 
(solid line with filled circles), and the exact result of the prolate spheroid (solid line). 
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